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Abstract 

The Mead- Fairfield geomagnetic field models have been 
used to trace field lines between the outer magnetosphere 
and the earth's surface. The results are presented in 
terms of ground latitude and local time contours projected 
to the equatorial plane and into the geomagnetic tail. Witl 
these contours various observations can be mapped along 
field lines between high and low altitudes. Low altitudes 
observations of the polar cap boundary , the polar cusp, 
the energetic electron trapping boundary and the sunward 
convection region are projected to the equatorial plane and 
compared with the results of the model and with each other. 
The results provide quantitative support to the earlier 
suggestions that the trapping boundary is associated with 
the last closed field line in the sunward hemisphere, the 
polar cusp is associated with the region of the last closed 
field line, and the polar cap projects to the geomagnetic 
tail and has a low latitude boundary corresponding to the 
last closed field line. The sunward convection region is 
associated with most of the magnetosphere outside the 
plasmapause , but whether the region extends all the way out- 
ward to. the last closed field line or inward to the 
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plasmapause is uncertain. Field lines from 6.6 R in 

E 

the geomagnetic equatorial plane intersect the earth 
between 67.4° and 65. 1°, with the exact value depending 
on the dipole tilt angle, season, local time and geo- 
magnetic activity. 
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Introduction 


Field and particle observations made along high 
latitude geomagnetic field lines form the basis for our 
current understanding of magnetosphere processes. Such 
measurements have generally been made either within 1 
of the earth's surface or at large distances near the 
equatorial plane. It is important to compare measurements 
in these two regions because the effects observed in one 
region may be controlled by physical processes in the oth.er 
region. Since most observed magnetosphere processes are 
related to charged particles which follow geomagnetic 
field lines, it is especially necessary to have accurate 
knowledge of the geomagnetic field in order to map the 
observed phenomena along field lines from one region to 
another. Within approximately 5 R £ of the earth a dipole 
model is usually adequate for such mapping purposes, but 
at greater distances , where external perturbations become 
larger relative to the internal field, models including 
these external sources must be used. 

A companion paper (Mead and Fairfield, 1973) presents 
four such field models which have been derived from over 
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4000 hours of vector magnetic field data obtained by four 
Explorer spacecraft. The present paper •will present results 
obtained from these models by tracing the field lines 
between the earth's surface and their intersection with either 
the equatorial plane or a cross section of the geomagnetic 
tail. Various statistical results obtained from low- 
altitude polar-orbiting spacecraft will be projected to the 
outer magnetosphere. This procedure will demonstrate use 
of the mapping technique and at the same time allow the 
earlier results to be viewed in a different context which 
may yield greater insight into physical processes. The 
low altitude measurements will be found to map to those 
regions of the magnetosphere that would be expected on the 
basis of various theoretical speculations. This fact can 
therefore be viewed as providing quantitative support for 

the theoretical expectations. Alternatively, one could 

( 

accept the validity of the theoretical ideas and argue that 
the mapping supports the accuracy of the new quantitative 
field models. Small differences between the projected 
experimental results and the expectations of various theories 
point out some critical outstanding questions of magnetospheric 
physics . 


4 



Field Line Tracing 


The work of Mead and Fairfield has produced four 

independent models corresponding to geomagnetic conditions 

that are quiet (K < 2) , disturbed (K 2) , super-quiet 

P P 

(K < o + ) and super-disturbed (K s 3) . Each of the models 
P p 

consists of a set of coefficients in a power series expansion 

for the disturbance field as a function of the position 

variables (solar magnetic coordinates) X,Y, and Z and tilt 

angle (the geomagnetic latitude of the sun) . The disturbance 

field when added to a dipole field yields a realistic 

magnetic field at any point in space and for any tilt angle. 

'j'hg range of validity for the models is restricted to 

distances < 17 R and the models have imposed north-south 

E 

symmetry and also east-west symmetry about a line 4° west of 

the sun, (ie. there is perfect north-south and east-west 

symmetry for zero tilts and for positive tilts the northern 

hemisphere has the same configuration as the southern 

hemisphere does for negative tilts) . It is important to 

realize that even though field models are available for 

four different K ranges, the models are static, whereas, 

P 

the actual magnetospheric configuration is dynamic and 
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often undergoes large changes on a time scale of tens of 
minutes. For this reason, the models are primarily of use 
in mapping phenomena where time variations have been 
averaged out. 

Since the magnetic field is available at all points 
of space, it is possible to trace field lines using any of 
the four models. Figure 1 presents some results obtained by 
tracing field lines and compares differences between the 
four models. At the top of Figure 1 is the equatorial cross- 
ing distance of the outermost closed field line in the noon 

meridian for zero tilt angles for each of the four models. 

The values range from 12.2 in the super-quiet model to 10.4 
in the super-disturbed model. It should be emphasized that 
these values emerge directly from the best fit models, in 
contrast to most earlier magnetospheric field models, where 
this distance is an adjustable parameter. The agreement 
with an observed magnetopause distance (Fairfield, 1971) of 
10.9 R is quite good, especially since the observed values 
include all tilt angles; the effect of non-zero tilts is 
to move the last closed field lines in the model closer to 
the earth. Figure 1 supports the common assumption that the 
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magnetopause corresponds to the last closed field line on 
the dayside of the earth. 

In the lower portion of Figure 1, the latitude where 
the last closed field line in the noon meridian plane inter- 
sects the earth is plotted for the four models. These 

o „ o 

values range from 80.5 at low K to 76.2 at high K , 

P P . 

which may be contrasted with the results of theoretically- 
derived field models (Williams and Mead, 1965; Choe and Beard, 
1973; see Roederer, 1969), which are generally greater than 
80° with realistic input parameters. The surface latitude 
of the last closed field line from the equatorial plane is 
independent of longitude in the solar hemisphere 

for any of the models. Thus the calculations support the 

concept (e.g. Aubry et al., 1970; Fairfield and Ness, 1970) 

that some lines of force that close on the day side of the 

earth under quiet conditions extend into the geomagnetic tail 

under disturbed conditions. The results are also consistent 

with those of Burch (1972) who found that the latitude of 

the cusp at low altitudes moves equatorward in association 

with a southward interplanetary field, which in turn is known 

to correlate well with K . 

P 
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Additional results from tracing three dimensional field 


lines are presented in Figure 2 (quiet model) and Figure 3 
(disturbed model) . These figures attempt to consolidate 
three dimensional information in a two dimensional figure 
by taking latitude and local time points on the earth's 
surface and tracing the corresponding field lines to the 
outer magnetosphere. The bottom portion of the figures 
show the intersection of these field lines with the 
equatorial plane, while the upper portions show their inter- 
section with a cross section of the tail at X = -10 R^. 

At 60°, where perturbations are small compared to the 
internal field, the latitude contours project approximately 
as circles, but at higher latitudes deviations from 
circularity increase. In Figure 2 the 69° latitude field 
line from midnight crosses the equatorial plane near X = -10 R 

ht 

with higher latitude field lines crossing the X = -10 R 

h» 

plane and going into the tail. For the disturbed model 
(Figure 3) the 69° field line goes further into the tail. 

The models are probably deficient at the larger Y values 
in Figures 2 and 3 at distances which exceed the 17 R„ 

hi 

limit of the original data. At these large distances the 
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Z component of the model field in the equatorial plane 
begins increasing in magnitude in a manner for which there 
is no observational support (Mead and Fairfield, 1973) . 

The X axis in Figures 2 and 3 are aligned along the earth -sutl 
direction and the east-west asymetry is due to the fact 
that the model is symetric about a direction 4° from the 
earth sun line. 

Mapping Low Altitude Measurements 

Figures 2 and 3 and their equivalents for the other 

2 models can be used for mapping, various low altitude 

measurements to the outer magnetosphere. Winningham (1972) 

used low altitude plasma measurements to determine the 

upper and lower boundaries of the polar cusp as a function 

of longitude for various K ranges. Evans and Stone (1972) 

P 

deduced the boundary of the polar cap by determing the low 
latitude termination to the polar cap access region for 

> 530 Kev solar electrons for two K ranges. If the 

P n . 

arguments are valid that (D the low energy polar cusp plasma 
has easy access to the magnetosphere on field lines that 
enter the magnetosheath and> (2) energetic electrons have 
access to open tail field lines in the downstream region, 
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these boundaries ought to be similar and compare favorably 

to the last closed field line determinations of the model. 

In Figure 4 the low altitude determinations of the polar 

cusp and the polar cap (the Evans and Stone measurements 

were averaged in two-hour local time intervals) during high 

K periods are projected to the outer magnetosphere along 
P 

the field lines of the super-disturbed model. The polar 

cusp which is observed in the 9:00 to 16:00 local time 

region is seen to project partially to closed field lines 

near the magnetopause and partially to the boundary of 

the geomagnetic tail (See also Akasofu et al , , 1973). The 

polar cap projects into the geomagnetic tail and actually 

overlaps the polar cusp to some extent, but no attempt has 

been made to overlay the two shadings. The low latitude 

boundary of the polar cap corresponds to the polar cusp 

region as well as could be expected considering the 

statistical uncertainties in the boundary determinations 

and the fact that the K ranges are somewhat different. 

P 

Both the cusp and the polar cap boundary data compare 
favorably with the predictions of the model. It is also 
interesting to note that the shape of the low latitude 
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polar cap boundary projected on the tail cross section 
corresponds very well to the high latitude boundary of 
the plasma sheet at 17 R_ (Bame et al., 1967). 

Figure 5 shows data similar to those of figure 4, 
only figure 5 is for low K conditions and the projection 

P -> 

is entirely into the equatorial plane. Again the polar 
cusp region is designated by shading but now the low latitude 
boundary of the polar cap is indicated by open circles. 

The additional data in figure 5 are those of McDiarmid 
and Burrows (1968). These authors were able to define 
two > 35 Kev electron boundaries in the auroral region 
and determine their average latitudinal position as a 
function of local time. The higher latitude boundary was 
defined as the position where flux levels reached background 
level . They suggested that this "background” boundary 
should be associated with the last closed field line 
of the magnetosphere. A lower latitude "smooth" boundary 
was defined as the location separating a low latitude 
region of stable fluxes from a higher latitude region 
of variable fluxes. They suggest that particles 
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equatorward of this boundary were completing longitudinal 
drift paths around the earth, whereas particles 
poleward of this boundary were undergoing pitch angle 
diffusion. The radial lines in Figure 5 represent 
equatorial plane projections of the model field lines 
corresponding to the boundaries of McDiarmid and Burrows. 

The solid lines represent field lines at the smooth 
boundary and the dashed lines those of the background 
boundary. The outer boundary in Figure 5 is the last 
closed field line of the quiet model and the inner oval 
is an average plasmapause of Chappell et al . , (1971). 

The agreement among the various types of data 
is quite good considering the uncertainties in the data. 

i 

However, the polar cusp is entirely on closed field 
lines in figure 5 and the last closed field lines of McDiarmid 
and Burrows lie somewhat inside of the last closed field lines 
of the model. It should be noted, however, that the model 
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field used corresponds to an average of 1 and if a 

slightly higher K model corresponding better to the particle 

P 

data had been available, more field lines would have gone 
into the tail (see Figure 1) and the agreement with the 
particle data would be improved. It should also be noted 
that the models have been used for zero tilt angle whereas 
the particle data comes from a variety of tilt angles. 

If non-zero tilt angles had been used in the model, 
fewer field lines would have closed on the day side and the 
agreement with the particle data would have been improved. 

It is interesting to note that there is a large dawn-dusk 
asymmetry in the trapping boundary of McDiarmid and Burrows, 
but only a very small asymmetry in the Evans and Stone data. 

A similar background field boundary of Fritz (1970) determined 
in a manner similar to that of McDiarmid and Burrows also 
did not show a large dawn-dusk asymmetry and it extended to 

somewhat higher latitude, presumably due to choice of a lower 

* 

background flux level. Asymmetries in the field model were 
precluded by imposing east-west symmetry on the model . Figure 
5 also implies that energetic solar electrons have access to 
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field lines which close as near the earth as 17 R . It 

E 

should be pointed out, however, that the model’s validity 
can be questioned in this region (Mead and Fairfield, 1973) 
and that a small change in the model could make this field 
line closing at 17 R extend many tens of earth radii into 
the tail , 

It should also be noted that Feldstein and Starkov 
(1970) suggest that the low latitude boundary of the auroral 
oval was associated with the > 35 Kev backgound boundary 
of McDiarmid and Burrows. This observation is consistent 
with the results of Frank and Gurnett (1971), which associate 
auroral arcs with "inverted V" particle events which occur 
at and above the > 45 Kev trapping, boundary. These results 
imply that aurora occur on "open" field lines, where the term 
"open" is defined in an operational sense in terms of the 
inability to support bouncing ~40 Kev electrons. The 
component of flux across the neutral sheet probably doesn’t 
go to zero near 20 R (e.g. Behannon, 1970) so it is likly 
that aurora actually occur on field lines that cross the 
neutral sheet at greater distances down the tail. 

Low altitude measurements of high latitude electric 
fields are a particularly interesting phenomena to project 
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to the equatorial plane since they are due to plasma motions 
in the outer magnetosphere. Heppner. (1972 a ) obtained 
measurements near the dawn-dusk meridian and Cauffman and 
Gurnett (1972) have presented data from a greater range of 
local times but from an experiment of somewhat lower 
sensitivity.* Although the measurements are of only one 
electric field component, both groups support a picture of 
sunward convection within the magnetosphere. When the 
effects of corotation are neglected and the low altitude 
measurements are projected to the equatorial plane they 
appear as in Figure 6. In this figure the two types of 
shading indicate the two data sets and the arrows the approx- 
imate direction of the convection. Also included is the 
projection of the Harang discontinuity (Harang , 1946) which 
has been suggested as the ionospheric locus along which occur 
the initial sudden magnetic field changes indicating the 
onset of expansion phases of magnetospheric substorms (Heppner, 
1972c) . 

The inner boundary of the convection region is slightly 
closer to the earth for the data, of Heppner, presumably due 
to greater sensitivity. Whether the convection region really 
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extends in as far as the plasmapause on a case-by-case basis 
is an important question which will have to be studied 
further . 

Another important question is whether the day side high 
latitude boundary of the convection region extends all the 
way to the magnetopause as proposed by Gurnett and Frank (1973) 
or whether anti-solar convection exists on closed field lines 
inside the magnetopause as is indicated in Figure 6. This 
important question will have to be resolved by direct 
measurements in the outer magnetosphere. Such anti-sunward 
convection inside the apparent magnetopause and near the 
equatorial plane has been observed under special circumstances 
on ATS-1 (Freeman et al . , 1968) and IMP-6 (Heppner, 1972b), 
but the generality of these observations has not been 
demonstrated. Recently Hones et al . (1972) and Akasofu et al . 

(1973) have observed anti-solar convection slightly to the 
nightside of the dawn-dusk meridian plane, with Vela satellites. 
This convection region is most prominent at higher latitudes 
but the lack of corresponding magnetic field measurements 
makes it particularly difficult to say if the convection 
occurs on open or closed field lines. Gurnett and Frank (1973) 
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also find that the high latitude boundary of sunward con- 


vection generally coincides with the high latitude termina- 
tion of > 45 Kev electrons. Comparing McDiarmid and Burrows 
statistical > 35 Kev electrons boundary (Figure 5) with the 
Cauffman and Gurnett data of Figure 6 does not reveal any 
large discrepancies which cannot be explained by statistical 
uncertainties in the data. 

Field Tracing at 6.6 R 

-L 

Figure 7 presents results obtained by tracing field 

lines from the synchronous satellite orbit at 6.6 R in the 

solar magnetic equatorial' plane. The data points indicate 

the geomagnetic latitudes where field lines originating 

at hourly intervals in the equatorial plane intersect the 

earth's surface. Three curves are shown for each model 

corresponding to a 30° tilt angle and tracing to the summer 

and winter hemispheres and a 0° tilt angle and tracing to 

either hemisphere. It can be seen that the surface 

latitude of field lines from the 6.6 R orbit exhibit a 

E 

variation with an amplitude which depends on the tilt angle 
and the season. The variation is largest in the summer for 
the most disturbed model (peak to peak amplitude of 2.3°) 
and smallest in the winter for the quietest model (peak to 
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peak amplitude of .3 ). It should be pointed out, however, 

that the curves of Figure 7 are not actual diurnal variations, 

since on any given day the tilt angle will vary whereas 

Figure 7 is plotted for a fixed tilt angle. The fact that 

the three points for each hour do not necessarily fall on 

a vertical line at the hour indicates that the longitude of 

the earth intersection also depends on the tilt angle. 

o 

The horizontal lines at 67.1 in Figure 7 represent the 
equivalent results for a dipole field. The maximum latitudes 
occur near noon where solar wind compression causes a field 
line from a higher latitude to intersect the 6.6 R orbit. 

hi 

At other longitudes the field is more streched out, causing 
lower latitude field lines to intersect this location. The 
fact that most of the points fall below the dipole value 
indicates that at most longitudes the ring current effects, 
causing inflated field lines, predominate over the compression 
effects of the solar wind. 

Summary and Discussion 

In summarizing the above results it can be said that 
there is general agreement between the predictions of the 
models and the particle measurements and their interpretation. 
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Specifically, the polar cap as determined from low altitude 
solar electron measurements projects to the geomagnetic tail 
and its low latitude boundary is close to the location of 
both the polar cusp and the limit of trapped > 35 Kev 
electrons. The polar cusp and the background trapping 
boundary correspond closely to the last closed field lines 
of the models. This fact is particularly significant since 
older field models generally give higher latitudes for last 
closed field lines and hence would not be in such good agree- 
ment with the particle clata. All phenomena have a similar 
variation with magnetic activity/ which are consistent with 
the transport of flux to the tail under disturbed conditions. 
The region of sunward convection as determined from low 
altitude measurements corresponds to almost the entire volume 
of the magnetosphere outside the plasmapause, but whether 
its boundaries extend all the way to the last closed field 
line or the plasmapause cannot be reliably determined 
from the available data. 

The general agreement between the prediction of the 
model and the results of particle measurements suggest that 
there are probably no grossly unrealistic features in the 
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field model and that there has been no serious misinterpretation 
of the particle data. More detailed comparisons between 
model predictions and magnetosphere measurements could, 
in principle, be used to answer many fundamental questions 
of magnetospheric physics. It would be interesting to know 
if polar cusp plasma resides entirely on field lines which 
enter the magnetosheath or whether the regions .designated 
by Winningham (1972) include particles that have diffused onto 
closed field lines, as has been suggested by Gurnett and 
Frank (1973). It is important to know if anti-solar 
convection exists on closed field lines inside the magneto- 
pause in order to understand physical processes at the 
magnetopause. Determination of whether or not sunward 
magnetosphere convection terminates at the plasmapause would 
indicate whether or not convection is important in the 
formation of the plasmapause. 

Although the search for the answer to these questions 
might seen straightforward, in practice the static nature 
of the models as opposed to the dynamic state of the actual 
magnetosphere dictate that meaningful comparisons be of a 
statistical nature. Unfortunately the number of particle 
measurements is limited, and therefore the results derived 
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from them have statistical uncertainties which preclude the 
detailed comparisons one might wish to make. One way to 
overcome this limitations would be to obtain additional 
particle measurements over more limited tilt and h.^ i anges 
and compare then with the appropriate model. A second method 
of answering some of these questions that have been raised 
by comparisons with field models is the direct comparison 
of simultaneous data from different experiments for various 
individual events. This approach has already been profitably 
exploited by Frank and Gurnett CL971 ) and should be continued 
and extended. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


- Equatorial crossing distance and earth inter- 
section geomagnetic latitude of the last closed 
field line in the noon meridian plane for dipole 
tilt angle of zero, plotted as a function of the 

median K value for each of the four geomagnetic 
P 

field models. 

- Contours of constant latitude (solid lines) and 

constant local time (dashed lines) at the earth’s 

surface as mapped along field lines of the 

MF73Q field model (K < 2) to the equatorial 

P 

plane (lower portion of figure) and to a tail 
cross section at X — -10 R (upper portion of 
figure) . The dipole tilt angle is zero. 

- Same as figure 2 but for the MF73D field model 

(K £ 2). 

P 





Figure 4 


Figure 5 


Figure 6 


- Low altitude experimental determination of the 

polar cusp and polar cap during periods of high 

K as mapped along field lines to the equatorial 
P 

plane (lower portion of figure) ^nd to a tail 
cross section (upper portion" of figure) . The 
super-disturbed model (K^ s 2) was used, with 
tilt = 0, The solid line on the dayside correspond 
to the last closed field line as determined by the 
model . 1 

- Low altitude measurements during quiet times as 
projected along field lines into the equatorial 
plane in the outer magnetosphere. The shaded 
region represents the polar cusp and the open 
circles represent the low latitude boundary 
of the polar cap. The dashed and solid lines 
emanating from the earth represent the equatorial 
plane projections of field lines corresponding 

to the 35 Kev background and smooth boundaries 
respectively . 

- Low latitude electric field measurements projected 
to the equatorial plane with the quiet model. The 
shaded areas indicate regions of sunward convection. 

*3 



Figure 7 


- Earth intersection geomagnetic latitudes of the 
field lines from various local times at the 
synchronous orbit in the geomagnetic equatoral 
plane. Four field models are shown for three 


tilt angles. 
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